Landau— Zener tunneling in the presence of weak intermolecular interactions in a 

crystal of Mn4 single-molecule magnets 
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A Mn4 single-molecule magnet (SMM), with a well isolated spin ground state of S = 9/2, is used 
as a model system to study Landau-Zener (LZ) tunneling in the presence of weak intermolecular 
dipolar and exchange interactions. The anisotropy constants D and B axe measured with minor 
hysteresis loops. A transverse field is used to tune the tunnel splitting over a large range. Using the 
LZ and inverse LZ method, it is shown that these interactions play an important role in the tunnel 
rates. Three regions are identified: (i) at small transverse fields, tunneling is dominated by single 
tunnel transitions; (ii) at intermediate transverse fields, the measured tunnel rates are governed 
by reshuffling of internal fields, (iii) at larger transverse fields, the magnetization reversal starts to 
be influenced by the direct relaxation process, and many-body tunnel events may occur. The hole 
digging method is used to study the next-nearest neighbor interactions. At small external fields, it is 
shown that magnetic ordering occurs which does not quench tunneling. An applied transverse field 
can increase the ordering rate. Spin-spin cross-relaxations, mediated by dipolar and weak exchange 
interactions, are proposed to explain additional quantum steps. 

PACS numbers: 75.50.Xx, 75.60.Jk, 75.75.+a, 75.45.+j 
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I. INTRODUCTION 

The nonadiabatic transition between the two states in 
a two-level system was first discussed by Landau, Zener, 
and Stiickelberg The original work by Zener con- 

centrated on the electronic states of a bi-atomic molecule, 
while Landau and Stiickelberg considered two atoms that 
undergo a scattering process. Their solution of the time- 
dependent Schrodinger equation of a two-level system 
could be applied to many physical systems, and it be- 
came an important tool for studying tunneling transi- 
tions. The Landau-Zener (LZ) model has also been ap- 
plied to spin tunneling in nanoparticles and molecular 

clusters 4,5,6,7,8,9,10,11,12^ 

Single-molecule magnetsi2*i^*i^ have been the most 
promising spin systems to date for observing quan- 
tum phenomena like Landau-Zener tunneling because 
they have a well-defined structure with well-characterized 
spin ground state and magnetic anisotropy^^-^^. These 
molecules can be assembled in ordered arrays where all 
molecules have the same orientation. Hence, macroscopic 
measurements can give direct access to single molecule 
properties. 

Since SMMs occur as assemblies in crystals, there is the 
possibility of a small electronic interaction of adjacent 
molecules. This leads to very small exchange interac- 
tions that depend strongly on the distance and the non- 
magnetic atoms in the exchange pathway. Up to now, 
such an intermolecular exchange interaction has been as- 
sumed to be negligibly small. However, our recent studies 
on several SMMs suggest that in most SMMs exchange 
interactions lead to a significant influence on the tun- 
nel processes. Recently, this intermolecular exchange in- 
teraction was used to couple antiferromagnetically two 
SMMs, each acting as a bias on its neighbor^^i^°i^^i^^i^'^. 



In this paper we present a detailed study of Landau- 
Zener tunneling in a Mn4 SMM with a well isolated spin 
ground state of 5* = 9/2. Using the standard and the 
inverse LZ method we show that spin-spin interactions 
are strong in SMMs with large tunnel splittings. By ap- 
plying transverse fields, we can tune the tunnel splittings 
from kHz to sub- GHz tunnel frequencies. We identify 
three regions depending on the applied transverse field. 
Next-nearest neighbor interactions, ordering, and spin- 
spin cross relaxations are studied. 

Several reasons led us to the choice of this SMM: (i) 
a half integer spin is very convenient to study different 
regions of tunnel splittings. At zero applied field, the 
Kramers degeneracy is only lifted by internal fields (dipo- 
lar, exchange, and nuclear spin interactions). A trans- 
verse field can then be used to tune the tunnel split- 
ting over a large range; (ii) Mn4 has a spin ground state 
S = 9/2 well separated from the first excited multiplet 
(S=7/2) by about 300 lO^; (ui) Mn4 has one of the 
largest uniaxial anisotropy constants, D, leading to well 
separated tunnel resonances; (iv) the spin ground state 
5' is rather small allowing easy studies of ground state 
tunneling; (v) Mn4 has a convenient crystal symmetry 
leading to needle shaped crystals with the easy axis of 
magnetization (being the c-axis) along the longest crys- 
tal direction; and (vi) the weak spin chain-like exchange 
and dipolar interactions of Mn4 are well controlled. 



II. STRUCTURE OF Mn4 AND MEASURING 
TECHNIQUE 

The studied SMM has the formula 
[Mn403(OSiMe3)(OAc)3(dbm)3], called briefly Mn4. 
The preparation. X-ray structure, and detailed physical 
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FIG. 1; (color online) (a) Unit cells of the Mn4 crystal. Only the cores of the Mn4 molecules are shown. The largest spheres 
are Mn, the smallest Si, and the others O atoms. The distances between next-nearest molecules are indicated, (b) Schematic 
drawing of the chain- like coupling between the Mn4 SMMs. (c) Scheme of the chain model where the S = 9/2 spin of each 
molecule is represented by an arrow. 



characterization are reported elsewhere-^'^. Mn4 crystal- 
lizes in a hexagonal space group with crystallographic 
C3 symmetry. The unit cell parameters are a — b ^ 
1.998 nm, c = 0.994 nm, a = /3 = 90°, and 7 = 120°. 
The unit cell volume is 3.438 nm'^ and two molecules are 
in a unit cell. The complex has a distorted cubane-like 
core geometry and is Mn3^^Mn^^. The C3 axis passes 
through the Mn^^ ion and the triply bridging siloxidc 
group (Fig. ^). DC and AC magnetic susceptibility 
measurements indicate a well isolated S — 9/2 ground 
stateS^. 

We found a fine structure of three in the zero-field res- 
onance fSect. HVB 1|) that is due to the strongest nearest 
neighbor interactions of about 0.036 T along the c— axis 
of the crystals. This coincides with the shortest Mn~Mn 
separations of 0.803 nm between two molecules along the 
c— axis, while the shortest Mn-Mn separations perpen- 
dicular to the c— axis are 1.69 nm and in diagonal direc- 
tion 1.08 nm (Fig. QJl). We cannot explain the value of 
0.036 T by taking into account only dipolar interactions, 
which should not be larger than about 0.01 T. We believe 
therefore that small exchange interactions are responsi- 
ble for the observed value. Indeed, the SMMs are held 
together by three H bonds C-H-0 which are probably 
responsible for the small exchange interactions. 

Fig. shows schematically the antiferromagnetic ex- 
change coupling between the Mn^^ {Si = 3/2) ions of one 
molecule and the Mn^^^ (6*2 = 2) ion of the neighboring 
molecule, going via three H bonds C-H-0 (not shown 
in Fig.njj). This leads to an effective ferromagnetic cou- 
pling between the collective spins {S = 9/2) of the SMMs 
(Fig. IB;) because the Mn^^ {Si = 3/2) ions and the Mn^^ 
{S2 = 2) ion in each molecule are antiferromagnetically 



coupled. 

All measurements were performed using an array of 
micro-SQUIDsS^. The high sensitivity allows us to study 
single crystals of SMMs of the order of 5 fim or larger. 
The field can be applied in any direction by separately 
driving three orthogonal coils. In the present study, the 
field was always aligned with the C3 axis of the molecule, 
that is the magnetic easy axis, with a precision better 
than 0.1°^-. The transverse fields were applied transverse 
to the C3 axis and along the a— axis. 



III. SPIN HAMILTONIAN AND 
LANDAU-ZENER TUNNELING 

The single spin model (giant spin model) is the sim- 
plest model describing the spin system of an isolated 
SMM. The spin Hamiltonian is 

n = -DSl - BSt + Htran. " Qf^Bf^oS ■ H (1) 

Sx, Sy, and Sz are the components of the spin operator, 
g « 2, /.tB the Bohr magneton; D and B the anisotropy 
constant defining an Ising type of anisotropy; Titrans, 
containing Sx or Sy spin operators, gives the transverse 
anisotropy which is small compared to DS^ in SMMs; 
and the last term describes the Zeeman energy associated 
with an applied field H. This Hamiltonian has an en- 
ergy level spectrum with (25" -t- 1) values which, to a first 
approximation, can be labeled by the quantum numbers 
TO = —S, —{S—1), S taking the z-axis as the quantiza- 
tion axis. The energy spectrum can be obtained by using 
standard diagonalization techniques (Fig.[2Jl. At = 0, 
the levels to = ±5" have the lowest energy. When a field 
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FIG. 2: (color online) Zeeman diagram of the 10 levels of 
the S = 9/2 manifold of Mn4 as a function of the field 
applied along the easy axis. The spin Hamiltonian param- 
eters are D = 0.608 K, B = 3.8 mK, and E = 0.032 K. 
The levels are approximately labeled with quantum numbers 
m = ±9/2, ±7/2, ...,±1/2. 



Hz is applied, the levels with m > decrease in energy, 
while those with to < increase. Therefore, energy levels 
of positive and negative quantum numbers cross at cer- 
tain values oiHz- Although Titrans produces tunneling, it 
can be neglected when determining the field positions of 
the level crossing because it is often much smaller than 
the axial terms. Without Titians and transverse fields, 
the Hamiltonian is diagonal and the field position of the 
crossing of level to with to' is given by 

_ n[D + B{m^ + m'^)] 

where n = —(to ± to') is the step index. 

When the spin Hamiltonian contains transverse terms 
(Titrans), the Icvcl crossiugs can be avoided level crossings. 
The spin S is in resonance between two states when the 
local longitudinal field is close to an avoided level cross- 
ing. The energy gap, the so-called tunnel splitting A, can 
be tuned by a transverse field (perpendicular to the 5*2 
direction) ^^V^i^^. 

The nonadiabatic tunneling probability Pm,m' between 
two states when sweeping the longitudinal field at a 
constant rate over an avoided energy level crossing was 
first discussed by Landau, Zener, and Stiickelberg iSi^. 
It is given by 

7rA2 , 

-Pm m' = 1 - exp "^^^ (3) 

2hgfiB\Tn — m'\fiodHz/dt 

Here, to and to' are the quantum numbers of the avoided 
level crossing, dHz/dt is the constant field sweeping rates, 
and h is Planck's constant. 

Fig. El presents two different methods to apply the LZ 
model: in Fig. the initial state is the lower energy 




magnetic field 



FIG. 3; (color online) Detail at a level crossing m with m' 
where the transverse terms (terms containing Sx or/and Sy 
spin operators) turn the crossing into an avoided level cross- 
ing. The initial state is the lower energy state in (a) (standard 
LZ method) whereas in (b) it is the higher energy state (in- 
verse LZ method). 

state (standard LZ method) whereas in Fig. it is the 
higher energy state (inverse LZ method). The tunneling 
probabilities are given by Eq.|3| In the simple LZ scheme, 
both methods should lead to the same result. However, 
when introducing interactions of the spin system with 
environmental degrees of freedom (phonons, dipolar and 
exchange interactions, nuclear spins, etc.), both methods 
are quite different because the final state in Fig. and 
the initial state in Fig. are unstable. The lifetimes 
of these states depend on the environmental couplings 
as well as the level mixing which can be tuned with an 
applied transverse field. We will see in Sect. IIII Al that 
comparison of the two methods allows the effect of envi- 
ronmental interactions to be observed. 



A. Landau— Zener tunneling in Mn4 

Landau-Zener tunneling can be seen in hysteresis loop 
measurements. Figs. 01a, and show typical hysteresis 
loops for a single crystals of Mn4 at several temperatures 
and field sweep rates. When the applied field is near an 
avoided level crossing, the magnetization relaxes faster, 
yielding steps separated by plateaus. As the temperature 
is lowered, there is a decrease in the transition rate due to 
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FIG. 4: (color online) (a) Hysteresis loop measurements of 
a single crystal of Mn4 at low temperatures (40 mK) where 
thermal activation to excited spin states can be neglected. 
The field is applied in the direction of the easy axis of magne- 
tization and swept at a constant rate between 0.002 and 0.14 
T/s. The dots labeled with Mpc are the magnetization after 
cooling the sample from 5 K down to 0.04 K in a constant 
applied field Hz- Mpc is used for the equilibrium magneti- 
zation Afoq in Sect. If V Al (b) Hysteresis loop measurements 
similar to Fig. 2 but at different temperatures and for a field 
sweep rate of 0.14 T/s. 



reduced thermally assisted tunneling. A similar behav- 
ior was observed in Mni2 acetate cluster3S2i21ii2ii22iS and 
other SMMa?'^i?^'^i?'?i?7i??. The hysteresis loops become 
temperature-independent below 0.4 K indicating ground 
state tunneling. The field between two resonances allows 
us to estimate the anisotropy constants D and B. We 
found: 
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FIG. 5: (color online) Minor hysteresis loops of a single crystal 
of Mn4. The magnetization was first saturated at -1.4 T. After 
ramping the field to zero at 0.14 T/s, the field was swept k- 
times back and forth (between 0.028 and 0.07 T in (a) and 
between -0.028 and 0.07 T in (b)) over the zero-field resonance 
with a sweep rate of 0.014 T/s. After the fc-th sweep, the field 
is quickly swept back to -1.4 T at a rate of 0.14 T/s leading to 
resonant tunneling at the level crossing (m, m') = ( — 7/2, 9/2) 
and (-5/2,9/2), and allow to determine ^oH^J"^ = 0.544 T 
and fioffi^^ — 1.054 T. The field interval of the k back and 
forth sweeps corresponds to zero reversed neighbor (0 RN, see 
Sect. If V B 11 1 in (a) whereas it goes over the RN, 1 RN, and 
2 RN transitions in (b). Note that the procedure in (a) leads 
to sharper steps and reduce spin-spin cross-relaxtion (SSCR) 
(Sect. IVJ because all reversed spins have two non-reversed 
neighbors. The transitions of SSCR are indicated in (b). 



where ffi^-* and iji^-* are the field positions of level cross- 
ings M = with S -1 and M = -S with S - 2. 

The influence of dipolar and intermolecular ex- 
change, which can shift slightly the resonance positions 
(Sect. IrVB can be avoided by performing minor hys- 
teresis loops involving only few percent of the molecules 
(Fig. [SJ. We found the field separations between the 
zero-field resonance and the first and second resonance 



are m ' = 0.544 T and m ' = 1.054 T. Using Egs. H 
andEl we find D = 0.608 and B = 3.8 mK. These values 
agree with those obtained from INS and EPR measure- 
ments^^. 

In order to explain the few minor steps (Fig. ^J, not 
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explained with the above Hamiltonion, spin-spin cross- 
relaxation between adjacent molecules has to be taken 
into account. Such relaxation processes, present in 
most SMMs, are well resolved for Mn4 because the spin 
is small (Sect-Ej). 

The spin-parity effect was established by measuring 
the tunnel splitting A as a function of transverse field 
because A is expected to be very sensitive to the spin- 
parity and the parity of the avoided level crossing. We 
showed elsewhere that the tunnel splitting increases grad- 
ually for an integer spin, whereas it increases rapidly for 
a half- integer spin^" . In order to apply quantitatively the 
LZ formula (Eq. |3Jl, we first checked the predicted field 
sweep rate dependence of the tunneling rate. The SMM 
crystal was placed in a high negative field to saturate the 
magnetization, the applied field was swept at a constant 
rate over one of the resonance transitions, and the frac- 
tion of molecules that reversed their spin was measured. 
The tunnel splitting A was calculated using Eq. 01 and 
was plotted in Fig. 2 of reference as a function of 
field sweep rate. The LZ method is applicable in the re- 
gion of high sweep rates where A_g/2,9/2 is independent 
of the field sweep rate. The deviations at lower sweep- 
ing rates are mainly due to reshuffling of internal fields^ 
(Sect. WJ^ as observed for the Fes SMM^. Such a be- 
havior has recently been simulated^iSi 



1. LZ tunneling in the limit of P^^^i << 1 

Fig. 3 of reference"^ presents the tunnel splittings ob- 
tained by the LZ method as a function of transverse field 
and shows that the tunnel splitting increases rapidly for a 
half-integer spin. Figs. 4a and 4b of reference present 
a simulation of the measured tunnel splittings. We found 
that either the second order term (£^(5*1 -|- 5^)) with E — 
0.032 K or a fourth order term {B^^[SX ^ 51)) with B44 
= 0.03 mK can equally well describe the experimental 
data. These results suggest that there is a small effect 
that breaks the C3 symmetry. This could be a small 
strain inside the SMM crystal induced by defects, which 
could result from a loss or disorder of solvent molecules. 
Recent inelastic neutron scattering measurements con- 
firm the presence of second and fourth order termgS. 



2. LZ tunneling for large probabilities Pm.m' 

Figs, ini and [7| show respectively the tunnel spfitting 
A and the LZ tunnel probabilities £±9/2 as a function 
of transverse field using the parameters of reference 
fSect. HTlAlll . P±9/2 increases rapidly to unity; for ex- 
ample, i-±9/2 — 1 for fftrans > 0.4 T and field sweep rates 
dH/dt smaller than 0.28 T/s. The Mn4 system is there- 
fore ideal to study different regimes of the tunnel prob- 
ability ranging from kHz to sub-GHz tunnel frequencies 
(Figs. EI). 




IQ-^-M 1 1 1 1 1 1—'- 
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(T) 

FIG. 6: (color online) Calculated tunnel splitting for Mn4 us- 
ing D = 0.608 K, B = 3.8 mK, and E = 0.038 K *^ The 
calculated A has been averaged over all possible orientations 
of the transverse field in order to represent the arbitrary ori- 
entation of the E term. The infiuence of nuclear spin broad- 
ening was taken into account by a Gaussian distribution of 
transverse field components with a half-width a — 0.035 T. 
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FIG. 7: (color online) Calculated LZ tunnel probabilities 
^±9/2 as a function of transverse field using the tunnel split- 
ting from Fig. |S| and the indicated field sweep rates. Only 
every 100-th calculated point is shown as a symbol. 

Fig. IHl presents the magnetization variation during LZ 
field sweeps for several transverse fields. The SMM crys- 
tal was first placed in a high negative field to saturate the 
magnetization and the applied field was then swept at a 
constant rate to the field value of -0.1 or 0.1 T for the 
LZ (Fig. ISk-d) or inverse LZ method (Fig. IHt-h) . At this 
field value, labeled 1, a transverse field was applied to in- 
crease the tunnel probability. Finally, the field is swept 
back and forth over the zero-field resonance transitions 
(m = ±9/2) and the fraction of molecules that reversed 
their spin was measured. 

Note severals points in Fig.|Sl (i) in Figs.|S^ and|Hl3, the 
magnetization increases gradually with each field sweep 
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FIG. 8: (color online) Magnetization versus applied field scans for (a-d) the LZ method and (e-h) the inverse LZ method. The 
indicated transverse fields i^trans are applied at point 1. 
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from M = -Ms to M = 0; (ii) in Fig.lHb, the field sweep 
from 2 to 3 shows first a decrease and then an increase 
of magnetization. This is due to next-nearest neighbor 
effects (Sect. HVBTJ; (iii) in Figs. |Ht and|Hli, the mag- 
netization increases (decreases) for a positive (negative) 
field scan. Note that the tunnel probability is 1 for trans- 
verse fields larger than 0.3 T (Figs. O, that is all spins 
should reverse for each field sweep; (iv) in Figs.|Ht and|Hfe, 
the magnetization increases much stronger for the field 
sweep from 1 to 2 than in Figs.|Hb and|Sl:; (v) in Fig.|Sli, 
the field sweep from 1 to 2 shows first an increase and 
then a decrease of magnetization; (vi) in Figs. IH^l to |Ht 
and Hi to ife, the magnetization tends to relax towards 
M = whereas in Figs.|Hli and|Hh, it relaxes towards the 
field cooled magnetization A/fc (Fig. EJi) . 

The result of a detailed study of the magnetization 
change AM for LZ field scans like those in Fig.|Slare sum- 
marized in Figs. 1^ and If 01 AM is obtained from AM = 
{Mf - Mi)^/\^\ where M; and Mf are the initial and 
final magnetization for a given LZ field sweep. Fig. [flU 
gives field sweep rate dependence for the field sweep from 
1 to 2. These graphs show clearly the crossover between 
the different regions presented in Fig. |H1 

We identify three regions: 

(i) at small transverse fields (0 to 0.2 T), that is 
P±9/2 « 1, tunneling is dominated by single tunnel 
transitions and AM follows the LZ formula (Eq.El). This 
regime is described in Sect. IIII A fl 

(ii) at intermediate transverse fields (0.2 to 0.7 T), that 
is tunnel probabilities P±q/2 between «0.f and wl, AM 
deviates strongly from Eq. |3| and is governed by reshuf- 
fling of internal fields; 

(iii) at larger transverse fields, the magnetization re- 
versal starts to be influenced by the direct relaxation 
process^** and many-body tunnel events may occur. 

The dominating reshuffling of internal fields in region 
(ii) can be seen when one compares AM in Figs. [SJa 
and ^ for the field sweep from 1 to 2 with those in 
Figs. ISt and |Hfe. A backward sweep gives a larger step 
than a forward sweep. This is expected for a weak ferro- 
magnetically coupled spin chain. Indeed, any spin that 
reverses shifts (shuffles) its neighboring spins to nega- 
tive fields. For a forward sweep this means that these 
spins will not come to resonance whereas in a backwards 
sweep, these spins might tunnel a little bit late during 
the field sweep. A more detailed discussion is presented 
in Sect. lfVBfl 

In region (iii) the direct relaxation procesa^^ between 
the two lowest levels starts to play a role. This can be 
seen by the fact that, during the application of the trans- 
verse field in point 1 (Figs.|Sli), the magnetization starts 
to relax rapidly. A direct relaxation process is indeed 
probable when the involved levels start to be mixed by 
the large transverse field. Because of this level mixing 
and the intermolecular interactions, multi-tunnel events 
are possible because neighboring spins start to be entan- 
gled. 

The inverse LZ method allows us to establish adiabatic 
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FIG. 9: (color online) The change of magnetization AAI = 
(Mf - Mi)^/\^\ as a function of transverse field fftrans 
for several LZ field sweeps; (a) LZ method and (b) inverse 
LZ method. Only every 30-th measured point is shown as a 
symbol. 

LZ transitions. Whereas for the standard LZ method 
the difference between an adiabatic and strongly deco- 
herent transition is difficult to distinguish, the inverse 
LZ method allows a clear separation. This is due to the 
fact that the equilibrium curve and the adiabatic curve 
are similar for the standard LZ method but not for the in- 
verse LZ method. For example. Fig. shows that there 
are more than 10 adiabatic LZ passages before the sys- 
tem reaches a disordered state. It is difficult to conclude 
this from Fig. ^ because strong decoherence would lead 
to a similar curve. 

It is important to note that the transition between re- 
gions (ii) and (iii) leads to the shoulder in Fig. 1 1 01 which 
should not be interpreted as quantum phase interference. 

Fig. 1^1 and Fig. IIOI are very rich with information and 
a complete understanding needs a multi-spin simulation. 

IV. INTERMOLECULAR DIPOLAR AND 
EXCHANGE INTERACTIONS 

SMMs can be arranged in a crystal with all molecules 
having the same orientation. Typical distances between 
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10: (color online) The change of magnetization AA/ — 
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FIG 

(Aff — Mi)^/\^ \ as a function of transverse field J/trans for 
the LZ field sweep from 1 to 2 and for the indicated field 
sweep rates; (a) LZ method and (b) inverse LZ method. 



molecules are between 1 and 2 nni. Therefore, inter- 
molecular dipole interactions cannot be neglected. An 
estimation of the dipolar energy can be found in the mean 
field approximation. 



dip 



An V 



(6) 



where V is the volume of the unit cell divided by the 
number of molecules per unit cell. Typical values of i^dip 
for SMM are between 0.03 and 0.2 K. More precise values, 
between 0.1 and 0.5 K, were calculated recently^^*^. 

In addition to dipolar interactions there is also the 
possibility of a small electronic interaction of adjacent 
molecules. This leads to very small exchange interac- 
tions that depend strongly on the distance and the non- 
magnetic atoms in the exchange pathway. Until recently, 
such intermolecular exchange interactions have been as- 
sumed to be negligibly small. However, our recent studies 
on several SMMs suggest that in most SMMs exchange 
interactions lead to a significant influence on the tunnel 
procesa^"i'"i'^^'i''\ 

The main difference between dipolar and exchange in- 
teractions are: (i) dipolar interactions are long range 
whereas exchange interactions are usually short range; 



(ii) exchange interactions can be much stronger than 
dipolar interactions; (iii) whereas the sign of a dipolar 
interaction can be determined easily, that of exchange 
depends strongly on electronic details and is very difficult 
to predict; and (iv) dipolar interactions depend strongly 
on the spin ground state S, whereas exchange interac- 
tions depend strongly on the single-ion spin states. For 
example, intermolecular dipolar interactions can be ne- 
glected for antiferromagnetic SMMs with S = 0, whereas 
intermolecular exchange interactions can still be impor- 
tant and act as a source of decoherence. 



A. Hole digging method to study intermolecular 
interactions 

Here, we focus on the low temperature and low field 
limits, where phonon-mediated relaxation is astronom- 
ically long and can be neglected. In this limit, the 
m — ±S spin states are coupled due to the tunnel split- 
ting A±s which is about 10""^ K for Mn4 (Sect. HhIi . 
In order to tunnel between these states, the longitudi- 
nal magnetic energy bias ^ = SMB'S'-ffiocai due to the lo- 
cal magnetic field iJiocai on a molecule must be smaller 
than A±s, implying a local field smaller than 10~^ T 
for Mn4 clusters. Since the typical intermolecular dipole 
fields for Mn4 are of the order of 0.01 T and the ex- 
change field between two adjacent molecules of the or- 
der of 0.03 T, it seems at first that almost all molecules 
should be blocked from tunneling by a very large en- 
ergy bias. Prokof 'ev and Stamp have proposed a solution 
to this dilemma by proposing that fast dynamic nuclear 
fluctuations broaden the resonance, and the gradual ad- 
justment of the internal fields in the sample caused by 
the tunneling brings other molecules into resonance and 
allows continuous relaxationii. 

Prokof'ev and Stamp showed that at a given longi- 
tudinal applied field Hz, the magnetization of a crystal 
of molecular clusters should relax at short times with a 
square-root time dependence which is due to a gradual 
modification of the dipole fields in the sample caused by 
the tunneling 



M(il„i) = A/in -f(Afeq(i?.) 



A/i 



OVW^ (7) 



Here M^n is the initial magnetization at time t = (af- 
ter a rapid field change), and M^qiH^) is the equilibrium 
magnetization at H^- Experimentally, Afgq is difficult 
to measure and we replaced it by the field cooled mag- 
netization Mfc{Hz) (Fig. 0^) Intermolecular exchange 
interactions are neglected in the theory of Prokof'ev and 
Stamp. 

The rate function TgqrtiHz) is proportional to the nor- 
malized distribution P{Hz) of molecules which are in res- 
onance at Hz 



sqrt 



Co 



A2 



Ed 4:h 



(8) 
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FIG. 11: Scheme of the hole digging method presenting the 
time dependence of temperature, applied field, and magneti- 
zation of the sample. 



where is the hne width coming from the nuclear spins, 
Ej) is the Gaussian half- width of P{Hz), and c is a con- 
stant of the order of unity which depends on the sample 
shape. Hence, the measurements of the short time relax- 
ation as a function of the applied field give directly 
the distribution P{Hz). 

Motivated by the Prokof 'ev-Stamp theorj*^, we devel- 
oped a new technique — which we called the hole digging 
method — that can be used to observe the time evolution 
of molecular states in crystals of nanomagnets^ and to 
establish resonant tunneling in systems where quantum 
steps are smeared out by small distributions of molecu- 
lar environmenl^. Here, it has allowed us to measure 
the statistical distribution of magnetic bias fields in the 
Mn4 system that arise from the weak dipole and exchange 
fields of the clusters. A hole can be "dug" into the distri- 
bution by depleting the available spins at a given applied 
field. Our method is based on the simple idea that after 
a rapid field change, the resulting short time relaxation 
of the magnetization is directly related to the number 
of molecules which are in resonance at the given applied 
field. Prokof 'ev and Stamp have suggested that the short 
time relaxation should follow a Vt— relaxation law [equa- 
tion lO]. However, the hole digging method should work 
with any short time relaxation law — for example, a power 
law 

M{H,,t) = Mi„ + (Afeq(i/.) - Afin)(r3hort(i?.)i)" (9) 



where Fshort is a characteristic short time relaxation rate 
that is directly related to the number of molecules which 
are in resonance at the applied field 7?z, and < a < 1 
in most cases, a = 0.5 in the Prokof 'ev-Stamp theory 
[equation jT))] and Fsqrt is directly proportional to P{Hz) 
(Eq. ISJ. The hole digging method can be divided into 
three steps (Fig. [TT|l : 



1. Preparing the initial state. A well-defined ini- 
tial magnetization state of the crystal of molecu- 
lar clusters can be achieved by rapidly cooling the 
sample from high down to low temperatures in a 
constant applied field H^. For zero applied field 
{Hz = 0) or rather large applied fields {H^ > 1 T), 
one yields the demagnetized or saturated magneti- 
zation state of the entire crystal, respectively. One 
can also quench the sample in a small field of a 
few milliteslas yielding any possible initial magne- 
tization Min. When the quench is fast (<< 1 s), 
the sample's magnetization does not have time to 
relax, either by thermal or by quantum transi- 
tions. This procedure yields a frozen thermal equi- 
librium distribution, whereas for slow cooling rates 
the molecule spin states in the crystal may tend 
to a partially ordered ground state. Sect. IIV B 21 
shows that, for our fastest cooling rates of 1 s, 
partial ordering occurs. However, we present a LZ- 
demagnetization method allowing us to reach a ran- 
domly disordered state. 

2. Modifying the initial state — hole digging. 

After preparing the initial state, a field -ffdig is ap- 
plied during a time idig, called digging field and 
digging time, respectively. During the digging time 
and depending on i?dig, a fraction of the molecu- 
lar spins tunnel (back and/or forth); that is, they 
reverse the direction of magnetization!— 



3. Probing the final state. Finally, a field HY°'°'^ is 
applied (Fig. Ill|l to measure the short time relax- 
ation Fshort (Eq- lU which is related to the number 
of spins that arc still free for tunneling after step 
(2). 



The entire procedure is then repeated many times but 
at other fields ifP™'^'' yielding Fshort (^^z, ^^dig, idig) which 
is related to the distribution of spins P(-ff2, -ffdig, idig) 
that are still free for tunneling after the hole digging. For 
^dig — 0, this method maps out the initial distribution. 

We applied the hole digging method to several samples 
of molecular clusters and quantum spin glasses. The most 
detailed study has been done on the Feg system. We 
found the predicted y/t relaxation (Eq. JTJ in experiments 
on fully saturated Feg crystals^SiS and on nonsaturated 
samples^. These results were in good agreement with 
simulationG^^i^^i^^'!^"^'^^t^. 
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FIG. 12: (color online) Typical square root of time relaxation 
curves for an Mn4 crystal measured at 40 mK. For each curve, 
the sample was first (a) saturated or (b) thermally annealed 
at = 0. Then the indicated field was applied and the short 
time relaxation of magnetization was measured. The slope of 
the lines gives Faqrt when plotted against the square-root of t 
as shown. 



B. Hole digging applied to Mn4 

Fig.Elshows typical relaxation curves plotted against 
the square-root of time. For initially saturated or ther- 
mally annealed magnetization, the short time square root 
law is well obeyed. A fit of the data to Eq. Q deter- 
mines Fsqi-t. We took Mcq = Mfc{Hz) of Fig. gji. A 
plot of Fgqrt versus H is shown in Fig. ^] for the sat- 
urated samples (Afin ~ M^), as well as for three other 
values of the initial magnetization which were obtained 
by quenching the sample from 5 K to 0.04 K in the pres- 
ence of a small field. The distribution for an initially 
saturated magnetization is clearly the most narrow re- 
flecting the high degree of order starting from this state. 
The distributions become broader as the initial magneti- 
zation becomes smaller reflecting the random fraction of 
reversed spins. However, a clear fine structure emerges 
with bumps at ±0.036 T and zero field which are due to 
flipped nearest neighbor spins fSect. HVB l|l . 



1 U 




Mo Hz (T) 

FIG. 13: (color online) Field dependence of the short time 
square-root relaxation rates Faqrt (J^z) for three different val- 
ues of the initial magnetization Min- According to equa- 
tion JSJ, the curves are proportional to the distribution P(Hz) 
of magnetic energy bias due to local internal field distribu- 
tions in the sample. Note the logarithmic scale for Fsqrt. The 
peaked distribution labeled Min ~ — 0.99Ma was obtained by 
saturating the sample, whereas the other distributions were 
obtained by thermal annealing. For Min << Ma, the curves 
are distorted by nearest neighbor effects. The peak at ±0.036 
and T are from molecules which have zero, one, or two 
nearest-neighbors (RN) molecules with reversed magnetiza- 
tion. 



Fig. I14r a') shows the short time relaxation rate for a 
digging field i?dig — 0.028 T and for several waiting times. 
Note the rapid depletion of molecular spin states around 
ffdig and how quickly the same fine-structure, observed 
in Fig. 1131 appears. The hole arises because only spins 
in resonance can tunnel. The hole is spread out because, 
as the sample relaxes, the internal fields in the sample 
change such that spins which were close to the resonance 
condition may actually be brought into resonance. The 
overall features are similar to experiments on a fully sat- 
urated Fes crystals'*^. However, the small chain-like in- 
termolecular interactions make the Mn4 system unique 
for a deeper study presented in the following. 



1. Chain-like intermolecular interactions 

The Mn4 molecules are arranged along the c— axis in 
a chain-like structure (Fig. ^ . The dipolar coupling be- 
tween molecules along the chain is significantly larger 
than between molecules in different chains. In addi- 
tion there are small exchange coupling between molecules 
along the chain (Sect.^Ji leading us to propose the fol- 
lowing model. 

Each arrow in Fig. ^] represents a molecule. The ± 
and — signs are the magnetic poles. The exchange cou- 
pling is represented by ± J. The ground state for a fer- 
romagnetic chain is when all spins are up or down with 
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FIG. 14: (color online) The field dependence of the short 
time square root relaxation rates rsqrt(-ffz) are presented on a 
logarithmic scale showing the depletion of the molecular spin 
states by quantum tunneling at i/dig = 0.028 T for various 
digging times tdig and A/in = —Ms- 

+ and — poles together and — J for all exchange cou- 
plings (Fig. 115b). For short, we say that all spins have 
zero reversed neighbors (2 RN). In order to reverse one 
spin at its zero-field resonance (to = —S and m' — S), a 
magnetic field has to be applied that compensates the in- 
teraction field from the neighbors. As soon as one spin is 
reversed (Fig. llSb ). the two neighboring spins see a pos- 
itive interaction field from one neighbor and a negative 
one from the other neighbor, that is we say for short that 
the two spins have one reversed neighbor (1 RN). The in- 
teraction field seen by those spins is compensated. Such 
a spin with 1 RN has a resonance at zero applied field and 
might reverse creating another spin with 1 RN (Fig. 115b ). 
The third case is when a spin has 2 RN (Fig. 11511). In this 
case, a negative field has to be applied to compensate the 
interaction field of the two neighbors. 

In summary, there are three possibilities for a given 
spin: RN, IRN, or 2RN with a zero-field resonance 
shifted to positive (0 RN), zero (1 RN), or negative (2 
RN) fields. The influence of the interaction fields of the 
neighboring molecules is taken into account by a bias 
field H"^^^^. The effective field acting on the molecule 
is therefore the sum of the applied field H^^^ and the 
bias field Hf'^'': 

1 ^ 

= iJ^PP -t- Hf^' = i/^PP + y J.sMk (10) 

where Mk is the quantum number of the neighboring 
molecule and Joff is an effective exchange coupling tak- 
ing into account of the nearest neighbor exchange and 
dipolar couphng. J^s ~ 0.01 K for Mn4. 

Because of the long range character of dipolar fields 
and the interchain dipolar couplings, the situation in a 
Mn4 crystal is more complicated. However, when the ex- 
change interaction is significantly larger than the dipolar 



FIG. 15: (color online) Schematical representation of a chain 
of spins with dipole and exchange coupling represented by -|- 
and — poles and ±J, respectively, (a) All spins have zero 
reversed neighbors (0 RN); in (b) and (c), the indicated spins 
have one reversed neighbor (1 RN), whereas in (d) it has 
two reversed neighbors (2 RN). The three cases of this two- 
neighbor-model lead to the fine structure of three for all quan- 
tum resonance steps. 



interaction, the long range character of the latter leads 
only to a broadening of the two-neighbor model. Fig. 1161 
presents schematically the distribution of internal fields 
of a randomly ordered, a partially ordered, and a com- 
pletely ordered state with zero total magnetization. Such 
distributions can be observed with a short-time relax- 
ation, presented in Figs.Elfor a Mn4 crystal with differ- 
ent initial magnetizations. 

We tested the two-neighbor model extensively using, 
for example, minor hysteresis loops and starting from an 
initially saturated state (Fig. [T3k ) . When sweeping the 
field over the zero-field resonance after a negative sat- 
uration field, resonant tunneling can only occur at the 
positive interaction field of 0.036 T. The corresponding 
step in M{H) indicates that few spins with RN re- 
versed creating spins with 1 RN. This leads to two steps 
when sweeping the field backwards over the zero-field res- 
onance, one for RN and one for 1 RN (Fig.lSja). When 
enough spins are reversed, a third step appears at 2 RN. 

Similar experiments can be done with the hole digging 
method (Sect.llV^. Digging a hole at the field of RN 
induces a peak at 1 RN (Fig. [Tl|l . 



2. Magnetic ordering in crystals of single-molecule magnets 

The question of magnetic ordering in molecular mag- 
nets has recently been addressed theoretically De- 
pending on the system ferromagnetic, antiferromagnetic, 
or spin glass like ground states with ordering tempera- 
tures between about 0.2 and 0.5 K have been predicted. 
Due to the slow relaxation of SMMs at low temperature, 
ordering might happen at non-accessible long time scales. 
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FIG. 16: (color online) Schema of the distributions of internal 
fields of (a) a randomly ordered, (b) a partially ordered, and 
(c) a completely ordered state with zero total magnetization. 
Here, A^^ and A''_ are the distributions for up and down spins, 
respectively. The fine structure with three bumps are due to 
three cases of zero reversed neighbors (0 RN), one reversed 
neighbor (1 RN), and two reversed neighbors (2 RN). 



Recent experimental studies concerned antiferromagnctic 
ordering in Feig SMM^, ferromagnetic ordering of high- 
spin molecules^, and partial ordering in the fast tunnel- 
ing regime of SMMsS^. We present here a simple method 
to show that partial ordering occurs in crystals of Mn4 
SMMs in the slowly tunneling regime. 

The first important step is to create a randomly disor- 
dered state fFig. 116b.). that is for any internal field value 
there are the same number of up and down spins. This 
means that for any applied field, no magnetization relax- 
ation can be observed because the tunneling from up to 
down is compensated by tunneling from down to up. 

We tried to achieve a randomly disordered state 
(Fig. 116b .) by a fast quench of the sample temperature 
from 5 K down to 0.04 K 1 s). When applying a small 
field, a -\/i-relaxation is observed (Fig. 112b ) showing the 
the sample was already partially ordered (Fig. 116( 3). 

We found that a randomly disordered state (Fig. ITCk ) 
can be achieved by sweeping back and forth the field over 
the zero-field resonance. During each sweep, few spins 
tunnel randomly back and fourth. When the Landau- 
Zener tunnel probability is small (Plz << 1), that is for 
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FIG. 17: (color online) Short time square root relaxation rate 
for an Mn4 crystal measured at 40 mK. (a) The ZFC curve 
was measured after thermal annealing at // = 0. For the 
other curves, a randomly disordered state with Min = was 
first created by sweeping back and forth the field over the 
zero-field resonance. After a waiting time tw at // = 0, the 
short time relaxation rate was measured. The field for one 
reversed neighbor (1 RN) is indicated. For longer tw, Fsqrt 
approaches the ZFC curve of a partially ordered state, (b) 
Similar to Fig. I17h . but during the waiting time a transverse 
field of 0.2 T was applied leading to faster ordering. 



fast sweep rates of 0.1 T/s for Mn4, and a large number 
of back and forth sweeps, a magnetization state can be 
prepared that shows only a very small relaxation when 
applying a small field (Fig. ITTk ) . Ordering can then be 
observed by simply waiting at = for a waiting time 
t„. The longer is tw, the larger is the relaxation, that is 
the distribution of internal fields evolves from a randomly 
disordered state (Fig. 116b ) to a partially ordered state 
(Fig. 116b ). In order to enhance the ordering, a transverse 
field can be applied during the waiting time (Fig. 117b). 
Note that we did not observe complete ordering (Fig. 116b) 
which is probably due to the entropy, similar to an infinite 
spin chain which will not order at T = due to entropy. 
It is also interesting to note that ordering does not quench 
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FIG. 18: (color online) Minor hysteresis loops of a single crys- 
tal of Mn4. The magnetization was first saturated at -1.4 T. 
After ramping the field to zero, the field was swept several 
times back and forth over a part of the tunneling step at the 
level crossing (m, m') = (—9/2,7/2): for curve (1) between 
0.38 and 0.62 T; for curve (2) between 0.38 and 0.5 T; for 
curve (3) between 0.52 and 0.62 T. After the back and forth 
sweeps, the field is swept to 1.4 T. The field sweep rate for all 
parts is 0.07 T/s. The fine structure due to nearest neighbor 
effects and SSCR are indicated. 



tunneling. 

V. SPIN-SPIN CROSS-RELAXATION IN 
SINGLE-MOLECULE MAGNETS 

We showed recently that the one-body tunnel picture 
of SMMs (Sect. is not always sufficient to explain 
the measured tunnel transitions. An improvement to the 
picture was proposed by including also two-body tun- 
nel transitions such as spin-spin cross-relaxation (SSCR) 
which are mediated by dipolar and weak exchange inter- 
actions between moleculesi^. At certain external fields, 
SSCRs lead to additional quantum resonances which 
show up in hysteresis loop measurements as well defined 
steps. A simple model was used to explain quantita- 
tively all observed transitionsiS Similar SSCR processes 
were also observed in the thermally activated regime of 
a LiYF4 single crystal doped with Ho ions^^ and for lan- 
thanide SMMs^o. 

In order to obtain an approximate understanding of 
SSCR, we considered a Hamiltonian describing two cou- 
pled SMMs which allowed as to explain quantitatively 13 
tunnel transitions. We checked also that all 13 transi- 
tions are sensitive to an applied transverse field, which 
always increases the tunnel rate. The parity of the level 
crossings was also established and in agreement with the 
two-spin modeli^. 

It is important to note that in reality a SMM is cou- 
pled to many other SMMs which in turn are coupled to 
many other SMMs. This represents a complicated many- 



body problem leading to quantum processes involving 
more than two SMMs. However, the more SMMs that 
are involved, the lower is the probability for occurrence. 
In the limit of small exchange couplings and transverse 
terms, we therefore consider only processes involving one 
or two SMMs. The mutual couplings between all SMMs 
should lead mainly to broadenings and small shifts of the 
observed quantum steps which can be studied with minor 
hysteresis loops. 

Fig. ^1 shows typical minor loops at the level cross- 
ing (m, m') = (—9/2,7/2). In curve (1), the field is 
swept forth and back over the entire resonance transi- 
tion. After about two forth and back sweeps, all spins 
are reversed. Note the next-nearest neighbor fine struc- 
ture that is in prefect agreement with the two-neighbor 
model fSect. irVB l|l . In curve (2), the field is swept forth 
and back over a SSCR transition (transition 7 in refer- 
ence — ) . Note that the relaxation rate is much slower 
because of the low probability of SSCRs and the fact that 
this transition is mainly possible for spins with RN or 
1 RN. In curve (3), the field is swept forth and back over 
a part of the level crossing (m,m') = (—9/2, 7/2) corre- 
sponding to spins with RN or 1 RN. In this case, the 
relaxation rate decreases strongly after the first forth and 
back sweep because the 2 RN spins cannot tunnel in this 
field interval. 



VI. CONCLUSION 

Resonance tunneling measurements on a new high 
symmetry Mn4 molecular nanomagnet show levels of de- 
tail not yet possible with other SMMs, as a result of 
higher symmetry and a well isolated spin ground state of 
S = 9/2. This has permitted an unprecedented level of 
analysis of the data to be accomplished, resulting in infor- 
mation not yet attainable with other SMMS. In particu- 
lar, Landau-Zener (LZ) tunneling in the presence of weak 
intermolecular dipolar and exchange interactions can be 
studied, using the LZ and inverse LZ method. The latter 
has not been applied to any other SMM. Three regions 
are identified: (i) at small transverse fields, tunneling is 
dominated by single tunnel transitions; (ii) at intermedi- 
ate transverse fields, the measured tunnel rates are gov- 
erned by reshuffiing of internal fields, (iii) at larger trans- 
verse fields, the magnetization reversal starts to be influ- 
enced by the direct relaxation process and many-body 
tunnel events might occur. The hole digging method is 
used to study the next-nearest neighbor interactions. At 
small external fields, it is shown that magnetic order- 
ing occurs which does not quench tunneling. An applied 
transverse field can increase the ordering rate. Spin- 
spin cross-relaxations, mediated by dipolar and weak ex- 
change interactions, are proposed to explain additional 
quantum steps. We would like to emphasize that the 
present study is mainly experimental, aiming to encour- 
age theorists to develop new tools to model the quantum 
behavior of weakly interacting quantum spin systems. 
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